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Abstract. 

The effect of an externally applied high frequency oscillating electric field on the 
critical nucleation field of superconductivity in the bulk as well as at the surface of 
a superconductor is investigated in details in this work. Starting from the linearized 
time dependent Ginzburg-Landau (TDLG) theory and using the variational principle 
we have shown the analogy between a quantum harmonic oscillator with that of the 
nucleation of superconductivity in bulk and a quantum double oscillator with that of 
the nucleation at the surface of a finite sample. The effective Hamiltonian approach of 
Cook et al [1] is employed to incorporate the effect of an externally applied highly 
oscillating electric field. The critical nucleation field ratio is also calculated from 
the ground state energy method. The results obtained from these two approximated 
theories agree very well with the exact results for the case of undriven system which 
establishes the validity of these two approximated theories. It is observed that the 
highly oscillating electric field actually increases the bulk critical nucleation field {Hc^ ) 
as well as the surface critical nucleation field (Hcs ) of superconductivity as compared 
to the case of absence of electric field (eo = 0). But the externally applied rapidly 
oscillating electric field accentuates the surface critical nucleation field more than the 
bulk critical nucleation field i.e. the increase of Hc-^ is 1.6592 times larger than that of 
Hr „ ■ 
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1. Introduction 

During recent years a lot of researcli activity is going on both in experimental and 
theoretical physics aiming at the understanding of dynamics of such systems which are 
exposed to strong time- dependent external fields [21 |3l HI [5] . Fundamental informations 
regarding high-temperature superconductors can be obtained from the high frequency 
electro dynamic response. Informations involving mixed state are extracted from these 
kind of studies P, [71 El El HDl E] ■ Also recent advances in microfabrication are creating 
new interesting opportunities for investigating the nucleation of superconductivity in 
type-II superconductor [I2l |13]. If one decreases the strength of an applied magnetic 
field below a certain critical value, a material can become superconducting and this 
critical value is known as nucleation field of superconductivity. Landau and Ginzburg 
[H] have shown that the value of this critical field for a bulk material {Hc^)^ equals k\/2 
{k is the dimensionless Ginzburg Landau parameter) times the value of its thermody- 
namical critical value (He). Saint James and de Gennes [15] discovered the existence of 
a higher critical field, H^^, by considering the nucleation at the surface of a semi-infinite 
material. Now the main question is, whether this ifcg is the universal upper limit of 
nucleation. In other words, is the ratio tt^ a universal constant ? In this perspective, 

-Hc2 

we investigate the high frequency nucleation field ratio jp- of type-II superconductor in 
the present paper. 

Time-dependent systems are generally more complicated than those of the corresponding 
time-independent ones. As a result, it is difficult to predict qualitative and quantitative 
behavior of driven systems even in cases in which it is very easy to understand the dy- 
namics of the corresponding time-independent ones. But there are certain methods by 
which these time-dependent systems can be described by an effective time-independent 
Hamiltonian [161 113 HI H [13 120] • This enables the qualitative as well as quantitative 
analysis of such driven systems more convincing. In order to consider the highly oscil- 
lating field we follow the sign convention of Denisov et al [2T] and the references therein. 
The Ginzburg-Landau theory for superconductivity represents one of the most useful 
tools available for the theoretical description of the nucleation of superconductivity in 
an applied field [22] . It starts with a free energy expansion, completely in line with the 
general Landau theory for condensed matter, with particular attention being paid on the 
gradient of the ordering quantity [23] . We use the linearized time dependent Ginzburg- 
Landau (TDGL) theory as the starting point of my discussion about the nucleation of 
superconductivity [SllESlEHlETlEHlEniEOlEIlES]. From the linearized TDGL theory 
we derive Schrodinger like equations as that of a single harmonic oscillator and a double 
oscillator for the bulk nucleation of superconductivity and surface nucleation of super- 
conductivity respectively. A. P. van Gelder have shown that the problem of nucleation 
resembles that of finding the ground state energy of a particle moving in a magnetic 
field and the ground state energy is inversely proportional to the nucleation field [33] . 
In the present study we want to demonstrate the link between the bulk nucleation field 
{Hc2) of superconductivity with that of finding the ground state energy of a single har- 



Nucleation of superconductivity under rapid cycling of electric field. 



3 



monic oscillator and the surface nucleation field (i^cs) of superconductivity with that of 
a double oscillator. We employ the effective Hamiltonian approach [16], [T71 [181 [D [13 [20] 
to incorporate the effect of highly oscillating field on the nucleation fields of supercon- 
ductivity. We calculate the nucleation field ratio tt^ through the ground state energy of 
a single harmonic oscillator and a double oscillator for the driven as well as non-driven 
cases. 

With this preceding background, we organize the rest of the paper as follows. In the next 
section, we discuss about the generalized linear TDGL theory of superconductivity. In 
this context we explore the connection between the bulk nucleation of superconductivity 
and the quantum harmonic oscillator and the similarity between a double oscillator and 
the surface nucleation of superconductivity. In section 3, we analyze the double oscil- 
lator in the presence of a high frequency field through the effective time-independent 
Hamiltonian method of Cook et al [Ij. By calculating the ground state energy of the 
driven single oscillator and the driven double oscillator, we determine the nucleation 
field ratio, -jp-, in the presence of a high frequency electric field. We summarize our 
findings and conclude in section 4. 

2. Linearized TDGL Theory & Nucleation of Superconductivity 

Long before the microscopic theory, a phenomenological approach to superconductivity 
was proposed by Ginzburg and Landau [22]. The idea was that the normal- 
superconducting transition is a thermodynamical second order transition. So one can 
apply to it the general theory of second-order transitions defining an order parameter 
in such a way that ip is zero in the disordered state (normal metal) and finite in the 
ordered state (superconducting metal). The free energy of a superconductor is given by 
[2l[32] 



F 

^ sn 



|2 , -^L/,|4 1 



alV'l' + ttI^I^ + ( - - '^^)^^]dv. (1) 



2 ' ' 2m* 

The transition from normal to superconducting state in a magnetic field is second order 
and near the transition point the order parameter is small, il) « ipoo and hence one 
can easily linerize the Ginzburg-Landau equation to the following form [211 1321 [3 



^(-zW-^l)V = -«^, (2) 

where ■0 stands for the complex superconducting order parameter and a is the first 
Ginzburg-Landau parameter, related to the temperature-dependent coherence length, 
^(T), by a = — 2^4^ (t) • '^^^ starting point of the theoretical description of nucleation 
{i.e. onset) of superconductivity in an applied magnetic field is this linearized Ginzburg- 
Landau equation (LGLE). One can easily identify that the Eq. (2) is identical with 
the Schrodinger equation for a free charged particle of mass m* and charge e* = 2e 
in a magnetic field H = V x A, with —a = \a\ playing the role of the energy 
eigenvalue. This property allows us to apply various familiar solutions and methods of 
usual quantum mechanics to the problem of nucleation in superconductivity. The lowest 
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eigenvalue of the LGLE gives the highest magnetic field at which the nucleation of the 
superconductivity can occur. Now the big question is what happens in nonstationary 
cases ? For instance, if one applies an external electric field which varies very slowly, is 
it possible that the order parameter will be given by the same equation as in the static 
case where the time is a parameter. On the other hand, if the field varies very rapidly 
with time will it be possible that the superconductor will respond to an average of the 
field as it happens in other systems in physics ? The last quary is our basic investigation 
in this work. Time dependent Ginzburg Landau (TDGL) model often gives a reasonable 
picture of superconducting dynamics [281 [29]. Unlike from its static counterpart, validity 
of the TDGL theory is much more limited. It is not enough just to be close to the critical 
temperature. The necessary condition is that the deviation from equilibrium is small; 
the quasiparticle excitations should remain essentially in equilibrium with the heat bath. 
It can normally be fulfilled for gap less superconductor [2ll|27j. So, we begin by writing 
down the time-dependent Ginzburg Landau equation that governs the dynamics of the 
superconducting order parameter [281 EH] : 

2 



^ dt fi ' 2m* 
Now choosing the electrical potential (p = — YXCos(a;t), H along the z axis with 
convenient gauge Ay = Hx and linearizing one can show 

— -^V +—^Hx—+[ — -^)x %l)-eeoxcos[ut)i) = ih— + \a\%l).{A) 
2m m*c oy ^ m*c"^ ^ -i ot 

Now using i/j = e "=^e^^^/(x)e~* and then following the method of Cook et al 

py one obtains 

+ - ^o)'/ = (|«| + — - TT-f /, (5) 



2m* dx"^ m*c^ ^ Am*uj'^ 2m 

1 hki, 

where xq = 



2.1. Bulk Superconductivity, Single Harmonic Oscillator and H^.^ 

In this subsection, we calculate the critical nucleation field of superconductivity in 
the presence of an external magnetic field for a large sample. We consider the large 
sample in the presence of a magnetic field H along the z axis and a convenient gauge 
is Ay = Hx. Equation (5) is our starting point of discussion about bulk nucleation. 
Equation (5) is same as that of a Schrodinger equation for a particle of mass m* bound 
in a harmonic oscillator potential with force constant The resulting harmonic 

oscillator eigenvalues are 

E^={n + l)huo={n+^-)h{^). (6) 
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In view of Eq. (5) these energy eigenvalues, En, are to be equated to (|«| + 4^^^ ~ ■ 
Thus 

H^-^^(H + ^-^). (7) 
(2n + l)e/i^ 4m*u;2 2m* ^ 

Here, we are concerned about the highest value of H i.e. Hc^ which is obviously given 
by the lowest eigenvalue (n = 0, £0 = and kz = Q). Thus 



rrold _ ^ C\a\ _ <P0 _ /Kjr /oN 



m c\a\ 

eft 2Tie 

where flux quantum 0o = |f • The above relation H°^^ = k^/2Hc gives us an important 
message. When > > He and vortex phase in the type-11 superconductor 

appears. On the other hand for < Hc^ < He and Meissner effect sets in, the 
mixed phase does not appear and one obtains type-I superconductor. 
On the other hand, for £0 7^ 0, 

2 2 

^ne. ^ rn*c{\a\ + £^) ^ ^ c ^ ^ ecel 

eh ehAuj^ AUuj^ 

It is evident from Eq. (9) that the highly oscillating electric field actually increases the 

2 

nucleation field by an amount in the bulk nucleation of superconductivity. 



2.2. Surface Superconductivity, Double Oscillator and H^ 



C3 



So far in my treatment of the Ginzburg-Landau equation at the mean-field level we have 
not taken into account the surface of the sample. At the surface of the superconductor, 
some additional boundary conditions need to be imposed on the solutions. One can quite 
reasonably expect that the presence of an interface between the superconductor and a 
non-superconducting material, such as a normal metal or an insulator must affect the 
nucleation of superconductivity in the material. We consider a specimen with a single 
plane surface and the external magnetic field to be parallel to the surface i.e. H = Hz. 
The superconducting sample is located in the half-space x > 0, while we take the non- 
superconducting material to be located in the half space a; < 0. The latter material 
is taken to be either vacuum or an insulating material. Then the superconducting 



/////////////////y/////////////////// 

► z 

H 



Figure 1. The finite sample in the upper half space a; > 0. 
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boundary condition imposed on ■0 in finite samples 



- ihV - —A 

c 



0. 



(10) 



It reduces to the Neumann boundary condition, Vi/'ln = 0, when the magnetic vector 
potential, A, can be chosen in a form with zero normal component at the boundary of 
the sample. In my case this becomes 

dil) 

x=0 



dx 



0, 



:iii 



where A = Hxy. We look for a solution of the form -0 = 
linearized TDGL equation (Eq. 4) with the constraints : 
df 

x=0 



dx 
df_ 

dx 







0. 



f{x) for the 

(12) 
(13) 



The complication arises because the boundary condition states that the solution must 
be flat at a position x = while the minimum is located at x = xq. When the minimum 
of the potential is located far from the surface {xq = oo) one can ignore the boundary 
condition, and when xq = the boundary condition is satisfied by the standard solution 
of usual Schrodinger equation of a harmonic oscillator. Thus in both the cases we 
obtain H = Hc^- One can easily understand that the surfaces have consequences for the 
solution to the LGLE only at the intermediate values of Xq, i.e. < Xq < oo- For the 
intermediate values of xq, we can think of solving the Schrodinger equation by employing 
the effective Hamiltonian method of Cook et al [1] and thus the effective Schrodinger 
like equation becomes 



-r d^f 2e^H^ 



{^x - Xoff 



2J2 



e e, 



-J = -a{xo)f. 



(14) 



2m* dx"^ m*c^ ^ 4m* cj^' 

This is an eigenvalue problem where the eigenvalue itself is xq dependent, and my task 
is to minimize this with respect to Xq subject to the boundary condition on / at the 
surface. Introducing x = xo = \f^^XQ. ^ = ^/^f^uj and (3 

can rewrite Eq. (14) as follows : 



2a 

huio ' 



one 



d'f 



+ (x-Xo)7-77^/ = /5/. 



(15) 



dx^ 4^]2 

Now my task is to find the lowest possible value of /3 subjected to the boundary 
conditions ^ = at x = (0,oo). This can be phrased as the following variational 
problem of minimizing the functional 



P 



lo^dx 



{x-Xo?P-'-^f 



(16) 



!^ dxP 

with respect to variations in /. The Euler-Lagrange equation for this variational problem 
is precisely the scaled differential Eq. (15). To do the minimization, we use the following 
trial wave-function : 

1, 



/(x) = exp 



(17) 
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With the help of this trial wave-function, we obtain 

r dx{{-bxr + ix' - 2XX0 + Xl) - SI} eM-bx') 



/o°°c?Xexp(-6x^) 



I 1 I XnA XO e'g.l 
4 ^ 46 ^ 2 y b 8^2 



2 + 26 + ^°~7?6"41]^- ^^^^ 



Now minimizing f3 with respect to xo, we obtain 
a/3 2 



+ 2X0 = (19) 
which yields xo = Substituting this back in /3 and again minimizing with respect 



dxo y/Trb 

to b, we find 



a/3 111 



and it gives us b — Jl — Substituting this back into P, we obtain 




Pmin-\ll---^. (21) 

From the definition of (3 one can relate 
Thus 

jjold jjold 

^ ~ /l _ 2 _ e2£2 ~ /-, _ 2 _ ece^ • (23) 
V n 4f^2 Y -L ^ 4?itj2H 

In view of equation (23) one can realize that the physically accessible region is defined 
as follows : Jl — ^ > > 0. Now solving equation (23) for 'H' one obtains 

TTold I ec£n Trnew 2 

Prom Eq. (24), we obtain for = 0, if^^f = 1.6592if°f and for So 0, the surface 
nucleation field Hg^"" = 1.6592//^/"'. The main message of Eq. (24) is that when the 
superconductivity occurs in a system, it starts to nucleate at the surface of an ideal 
defect-free sample and not in the interior of the sample. If the sample has defects in the 
interior, superconductivity starts to nucleate in the vicinity of such defects. On the other 
hand the high frequency field further accentuates the nucleation at the surface rather 
than in the interior of the sample which is evident from equation (24). For the bulk 

2 

sample the increasing amount due to the rapidly oscillating field is (see equation 
9). On the other hand the rapidly oscillating field increases the surface nucleation field 
by an amount 1.6592^^ (see equation 24). Thus the rapidly oscillating electric field 
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accentuates the surface nucleation field by 1.6592 times more than the bulk nucleasion 
field increasing factor. From equation (9) and equation (24) one can determine the real 
values for this enhancement in nucleation field. If an experimentalist uses Eq = 0.01 
V/m, u = 10 GHz, the enhancement of the bulk nucleation field is 0.114 T more than 
the without driven result. For the same values of electric field strength and frequency, 
the enhancement of surface nucleation field is 0.189 T more than the non-driven system. 
Before concluding this subsection we want to show the analogy between the problem of 
nucleation of superconductivity at the surface with that of a double oscillator. Equation 
(14) is still have the form of Schrddinger equation for a particle in a harmonic well cen- 
tered at xq, but the boundary condition, Eqs. (12) and (13), means that the eigenvalue 
depends crucially on the value of Xq. One can incorporate the boundary condition by an 




Figure 2. (color online) (a) Surface and interior nucleation at with single oscillator 
potential form; (b) Surface nucleation at Hc-^ with double-oscillator potential form. 



image method. We consider a particle moving in the potential V{x) of Fig. 2(b) which 
is my original harmonic well together with its reflection at the surface. The ground-state 
wave function iPq{x) in V{x) must be symmetric about x = 0, as required in Eq. (12), 
and for x > it satisfies Eq. (14). Thus ipo{x) for x > is the solution to our problem, 
and the corresponding eigenvalue Eq gives the critical field. We can compare Eq for 
various xo with the eigenvalue E in the single harmonic well. For xq ^ oo, Eq ^ E, 
and for xq = 0, Eq = E again. On the other hand for intermediate values of xq, Eq is less 
than E because V{x) is smaller than the single harmonic potential in some region. The 
new surface eigenfunction must have a lower eigenvalue than the interior ones because 
it arises from a potential that is lower and broader than the single harmonic well about 
Xq. Thus one can say that the problem of nucleation of superconductivity at the surface 
is analogous to that of double oscillator problem. 
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3. Driven Double and Single Oscillator : Critical Field Ratio 

The analysis of the previous section has clearly demonstrated that the problem of 
nucleation of superconductivity in the bulk and at the surface are analogous to that 
of a single harmonic oscillator and double oscillator respectively. A. P. van Gelder had 
shown that the problem of nucleation resembles that of finding the ground state energy 
of a quantum particle and the nucleation field is inversely proportional to the ground 
state energy [33]. Henceforth, one can easily calculate the nucleation field ratio, 
by calculating the ground state energies of a single harmonic oscillator and a double 
oscillator. Due to the recent developments in the regime of superconducting high 
frequency devices [351 ES] , we would like to exhibit the effect of high frequency field on the 
nucleation of superconductivity in this section. The effect of high frequency field can be 
taken care of by simply following the effective Hamiltonian method [T6 l fTT l flS l fTl [T91I20] . 
By following Cook et al [Ij, we convert the time-dependent problem into an effective 
time independent one and then calculate the ground state energy by following standard 
stationary quantum mechanics procedures. Thus the main objective of this section is to 
find the nucleation field ratio, tt^, from the eigenvalue solutions of the single harmonic 
oscillator and the double oscillator in both the cases i.e. in the presence and absence of 
the high frequency field. 

We consider a quantum particle moving in a double oscillator potential Vo(|x|) = 
|m*ct;Q(^|x| — Xq) and driven by a high frequency monochromatic force. Thus the driven 
double-oscillator Hamiltonian becomes 

Ti = — 1 — m*a;Qf|x| — Xq) — axcos(ut). (25) 

2m* 2 ^ ^ 

Now following Cook et al [1\ we can express the effective time-independent potential 

1 2 Ci^ 

Veffix) = -m*ul[\x\ - xo) - (26) 

To determine ground state energy of this driven double oscillator, we consider the 
normalized ground-state eigenf unctions of single oscillators centered at ±Xo : 

«^ = (^)^e-(^^^«)', (27) 



h 



where we have introduced dimensionless variables j = rjx and 70 = r]XQ with rj 
Now, we compute the expectation value of the Hamiltonian (7Y), in the state -00 
-^(u+ + U-). The norm of ipo is A^o = {ipo\'4^o) = 1 + e~'''o. To compute the expectation 

2m* dx'^ ■ 



value of the kinetic energy operator (T = —-^-f^) we need to use the identity 



2 

Tm± = ^K-(7T7o)V]. (2^ 



Thus, we find 

(t/^oI^o) 4m* V eTo + 1 
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To obtain the expectation value of the potential energy operator V = ^m*ujQ(^x^ — 
2xo\x\ + Xg) — ^^r^2 ; we need to compute {ijjo\x'^\ipQ) and (^/'ol k| |'0o)- So, we obtain 



Finally the ground state energy of the driven oscillator is given by 



(30) 




In Fig. 3, we plot the ground-state energy versus 70 for the double-oscillator in the 
presence as well as absence of the external high frequency field. As pointed out by A. 
P. van Gelder the critical nucleation field is inversely proportional to the ground state 
energy. Thus, we calculate the critical nucleation field ratio (jf^) from the numerical 
computation of the minimum of ground state energies of the double and single oscillators 
(Eq. (31)). Minimum ground state energies for the double oscillator (dose) in the 
presence and absence of the external field are computed from Eq. (31) by putting 
a = 1.0, uj = 4.0, uo = 0.5 and a = respectively. On the other hand the ground state 
energies for the single oscillator (sosc.) in the presence and absence of the field (a = 0) 
can be computed from the Eq. (31) by simply putting 70 = 0. From our computation, 
we obtain the minimum of ground state energies for the double-oscillator with external 
rapidly oscillating force (wf) and without force (wof) cases are (-ff)o,mm = 0.2808 and 
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{H)q'^^^ = 0.2964 respectively. Thus the ratio of the critical fields for the force-free case 
is given by 

^ = Sflh^ = = 1.6869. (32) 

^cf {H)'ZL. 0-2964 

We have to mention that a more careful numerical analysis based on the hypergeometric 
functions in the non-driven system yields (exact result) 

^old 

WTd = 1-695. (33) 

From the variational principle one obtains H"lf' = 1.6592i7°2'^ and from the ground state 
energy method, we obtain Hf^'^ = 1.6869i/°f . Now comparing these two results with 
equation (33) one can say that the results obtained from the two approximated theories 
(variational principle and ground state enrgy) agree very well with the exact results. 
This establishes the validity of these two approximated theories. 

Now, we employ the ground state energy method to deternine the nucleation field ratio 
in the presence of the externally applied high frequency field : 

^^(^^ 04687 
{H)o,Lc. ('■2652 

In view of Eq. (34) and Eq. (32), we can conclude that the nucleation field ratio for the 
driven system is slightly larger than the non-driven system. The basic inference that 
one can draw from the above discussion is that the high frequency field accentuates the 
surface critical nucleation field of superconductivity more than that of the bulk critical 
nucleation field. 



4. Summary & Conclusions 

In this section, we briefiy summarize all my derived results and then conclude this 
paper. We investigate the effect of high frequency electric field on the nucleation of 
superconductivity in the interior of a large sample as well as at the surface of a finite 
sample. Employing the linearized time dependent Ginzburg-Landau theory for this pur- 
pose we have shown the analogy between quantum single harmonic oscillator and the 
nucleation in the interior of a large sample. The similarities between the nucleation of 
superconductivity at the surface and a quantum double oscillator is also demonstrated. 
We invoke two approximate theories to derive critical nucleation field ratio jf^. The 
first approximate theory is based on the variational principle and is solved analytically. 
On the other hand the second approximate theory is based on the ground state energy 
method of A. P. van Gelder [33]. We determine the minimum of the ground state energy 
by the numerical method. Since critical nucleation field is inversely proportional to the 
ground state energy one can easily compute the critical nucleation field ratio. The va- 
lidity of these two approximated theories is checked by comparing the results obtained 
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from these two theories with that of the exact results for the case of nondriven system 
and they agree very well with each other. The effect of high frequency oscillating time- 
periodic field is taken care of by the effctive time-independent Hamiltonian method of 
Cook et al [V\. It is observed that for electric field strength e = 0.01 V/m and frequency 
a; = 10 GHz one can obtain 0.114 T and 0.189 T much more enhancement in the bulk 
nucleation and surface nucleation field respectively. 

In conclusion, we examine in details the effect of high frequency field on the nucleation 
of superconductivity. Using the variational method and ground state energy method 
we have shown that the high frequency field actually accentuates the surface critical 
nucleation field of superconductivity more than the bulk critical nucleation field. The 
enhancement of the surface critical nucleation field is 1.6592 times more than the en- 
hancement of the bulk critical nucleation field. One can now say that the ratio jf^ is 
not universal i.e. i/c3 is not the universal upper limit of nucleation. One can obtain 
higher than Hc^ field by applying high frequency oscillating electric field. Our results 
will be helpful in analyzing recently developed high-frequency superconducting devices 
[351 [36]. 
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